arXiv: 1505.0457lv4 [hep-ph] 31 Mar 2016 


Implication of the observed ^ ppn^ for studying the pp i^(3n0)n^^ process 


Hao Xu'’^* Ju-Jun and Xiang Liu'’-4« 

^Research Center for Hadron and CSR Physics, Lanzhou University & Institute of Modern Physics of CAS, Lanzhou 730000, China 
^School of Physical Science and Technology, Lanzhou University, Lanzhou 730000, China 
^Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 
State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, Chinese Academy of Sciences, Beijing 100190, China 

We study the charmonium pp ifi(3n0)7t° reaction using effective lagrangian approach where the contribu¬ 
tions from well established N" states are considered, and all parameters are fixed in the process of e^e^ —> ppr/’ 
at center of mass energy ^fs = 3.773 GeV. The experimental data on the line shape of the mass distribution of 
the e*e^ —> ppn^ can be well reproduced. Based on the studying of e*e^ —> pptt^, the total and differential cross 
sections of the pp ^/^(3770)7r° reaction are predicted. At the same time we evaluated also the cross sections of 
the pp —> i/'(3686);r® reaction. It is shown that the contribution of nucleon pole to this reaction is largest close 
to the reaction threshold. However, the interference between nucleon pole and the other nucleon resonance can 
still change the angle distributions significantly. Those theoretical results may be test by the future experiments 
at PANDA. 

PACS numbers: 13.20.Gd, 13.75.Lb 


I. INTRODUCTION 

As a forthcoming facility in future, the Anti-Proton Anni¬ 
hilations at Darmstadt (PANDA) experiment will focus on the 
production of charmonium, which is govern by nonperturba- 
tive effect of quantum chromodynamics (QCD) [1]. Before 
PANDA run, there were pioneering theoretical studies of the 
charmonium production in the pp annihilation processes [2- 
9]. By calculating two hadron-level diagrams introduced by 
the Born approximation, Gaillard and Maiani firstly studied 
the differential cross section of the charmonium production 
plus a soft pion in the pp reaction [2]. In Ref. [3], the cross 
sections of the chamonium (*10 production accompanied by a 
light meson (m) from the process of pp ^ 'V + m was cal¬ 
culated by combing with the measured partial decay widths 
of charmonium decay into ppm. And then, Barnes and Li 
proposed an initial state light meson emission model for the 
near threshold associated charmonium production processes 
PP —> ('P = rjc, XcO, Xci), and the total and dif¬ 

ferential cross sections for these reactions were evaluated [4— 
6]. It is also found that the cross section of pp —> 7r°'P near 
threshold may be affected by the Pauli 7/iApp coupling [5]. 
Furthermore, Lin, Xu and Liu revisited the issue of the pro¬ 
duction of charmonium plus a light meson at PANDA, where 
the contribution of form factors (FFs) to these processes are 
included [7]. Recently, Pire et al. studied the associated 
production of a 7/t/r and a pion in antiproton-nucleon anni¬ 
hilation in the framework of QCD collinear factorization [8], 
while in Ref. [9], the exclusive charmonium production pro¬ 
cess pp ^ TT^Jftftwas studied within a nucleon-pole exchange 
model by including off-shell hadronic FFs and a complete 
Lorentz structure with a ppJl^i Pauli strong coupling. The 
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contributions from the intermediate N* states are also studied 
in Ref. [9], and it was found that one can not ignore the con¬ 
tributions of the N* resonances in the pp ^ 7r‘’7/i/t reaction. 

The experimental activity on the charmonium decays have 
run in parallel. These decays are of interest because they can 
be used to study the associated charomonium production in 
pp annihilation. In 2014, the BESIII Collaboration reported 
the analysis of e^e" in the vicinity of i/^(3770) [10]. 

In addition to the Born cross section of e^e~ PP^^, the 
corresponding pn^ and pit^ invariant mass distributions of 
PpTt^ process are also measured [10]. These new ex¬ 
perimental information in Ref. [10] allows us to further per¬ 
form a comprehensive study of —> if/{377Q) —> pptfi, 

which stimulates our interest to study the contribution of ex¬ 
cited nucleon resonances (N*) to e^e~ i/^(3770) ^ ppn^ 
and i/^(3770) production from pp ^ ^{377ff)n^ reaction. 

The nucleon is the simplest system in which the three col¬ 
ors of QCD can combine to form a colorless object, thus it is 
important to understand the internal quark-gluon structure of 
the nucleon and its excited N* states, and the study of excited 
N* states is an interested research field of hadron physics [11], 
which can make our knowledge of hadron spectrum abundant. 
A very important source of information for the nucleon inter¬ 
nal structure is the N* mass spectrum as well as its various 
production and decay rates, while the charmonium decay into 
ppn^ is an ideal platform to study excited N* nucleon reso¬ 
nances, because it provides an effective isospin 1/2 filter for 
the nN system due to isospin conservation [12-14]. 

In this work, we introduce excited N* nucleon resonances 
in the process of —> ift(377Q) — > pptfi. By fitting the 

p7r° and p7r° invariant mass distributions of the cross section 
of pptt^\ we extract the information of couplings of 

N*N7t and ij/(3770)N*N, which not only reflects the inner fea¬ 
tures of discussed N*, but also helps us to learn the role played 
by N* in the i/'(3770) ^ pptt^- 

Based on our studying on the iJ/{377Q) —> ppTr*’ 

process, we move forward to study the pp ^ jA(3770)7r° 
reaction, which is due to the cross relation between the 
1/^(3770) ^ ppn^ decay and the pp ^ ijj{3770)n^ reac- 
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tion [13]. Here, these extracted parameters from our study 
of e"e~ i/^(3770) —» ppn^ will be employed to estimate 
the production rate of pp i//{371Q)n^ and relevant features. 
We calculate the total and differential cross sections of the 
pp iI/{311Q)7jP reaction. It is shown that the contribution 
of nucleon pole to this reaction is the largest close to the re¬ 
action threshold. However, the interference between nucleon 
pole and the other nucleon resonance affects significantly and 
could change the angle distributions clearly. Additionally, 
there were abundant experimental data of i/^(3686) —» ppn^ 
given by BESIII [14], where BESIII released the branching 
ratio B(iA(3686) ^ ppn'-^) = (1.65 + 0.03 + 0.15) x lO^^ 
and the measured pn^ and pif^ invariant mass spectra [14]. 
This experimental status related to i/r(3686) makes us extend 
the above study to the i/r(3686) ^ ppn^^ decay, and also the 
pp i/'(3686)7r‘^ reaction. Our studies provide valuable in¬ 
formation to future experimental exploration of i/r(3770) and 
1/^(3686) productions plus a pion through the pp interaction at 
PANDA. 

This paper is organized as follows. After introduction in 
Sec. I, we present the detailed study of PP^^ by in¬ 

cluding the excited N* nucleon resonances (see Sec. II). In 
Sec. Ill, we further calculate pp i/'(3770)7r“ by combining 
with these results obtained in Sec. II. In Sec. IV, we adopt 
the similar approach to study i/f(3686) —» ppn'^ decay and the 
pp ^{36%6)n^ process. The paper ends with a discussion 
and conclusion. 


II. EXCITED A* NUCLEON RESONANCE 
CONTRIBUTIONS TO ^ 1^(3770) ^ 

Eirst, we study the process —» ppn^ with an ef¬ 

fective Lagrangian approach. In hadron level, the process 
e^e~ PP^^ in the vicinity of ^(3770) is described by the di¬ 
agrams shown in Eig. 1 . In Eig. 1 (a), and annihilate into 
photon, which couples with charmonium t//{377Q). And then, 
i/^(3770) interacts with final states, where we consider the con¬ 
tributions from nucleon-pole (s Pn) with 7^ = ^ and five 
N* states that are well established [15]: At(1440) (s Pn) with 
jP = i^, A(1520) (= Di 3 ) with JP = A(1535) (= Sn) 

with jP = A(1650) (= 5ii) with and A(1720) 

(s P 13 ) with - 1^. Additionally, we also consider the 
background contribution, where the annihilation directly 
into ppn without intermediate i/f(3770), which is shown in 
Fig. 1 (b). 

To compute the contributions of these terms, we use the 
effective interaction Lagrangian densities for each vertex. Eor 
the y^(377Q) coupling, we adopt the vector meson dominant 
(VMD) model, where a vector meson couples to a photon is 
described by [16] 
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FIG. 1: (color online). The Feynman diagrams for the process 
—> ppTT^ in the vicinity of ^('(3770). 


e!fv can be fitted through V e^e : 
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where |/^ = (M^ - 4 ot2)F2/2 ^ My 12 is three momentum 
of an electron in the rest frame of the vector meson, a - 
e^KAji) = 1/137. Using B(^(3770) ^ e+e") = (9.6 + 0.7) x 
10“^ [15], we obtain ejf^ono) = 0.0053. 

The JftfrNN and NNn couplings are described by: 


£kNN = -TT —■ d^nN, (3) 

■C^NN - -gif/NNNyijVPN, (4) 


where stands for the vector field of i/r(3770). We take 
Sknn - 13.45. 

For the N*Nti and ^N*N vertexes, we adopt the Lagrangian 
densities as used in Refs. [17-23]: 


-CjiNPn 

Nysy^T ■ dPjiRp^^ + h.c., 

IfflN 

(5) 

■CnNSu 

- -gnNSn^T ■ TlRsn + h-C., 

( 6 ) 

■LjiNPu 

_ _ SkNPu ^ duTiR'i + h.c., 

niN ^ 

(7) 


= 2 Nysy^T ■ d^dynRl, + h.c., 

( 8 ) 




-C.ipNP\\ 

= -g^NPnNypV^Rp,, + h.c.. 

(9) 

J^ijjNS 11 

= -gtpNSn^yiypy^Rsu+h-C., 

( 10 ) 


= -ig,pNPuNysV^R^p^^+h.c., 

( 11 ) 


= -gfNDijNVpR^P,^^ + h.c.. 

( 12 ) 


where P is a A* field. 

For the intermediate nucleon-pole or N* state, a Breit- 
Wigner form of its propagator Gj(q) can be written as [24] 


In above expression. My and fy are the mass and the decay 
constant of the vector meson, respectively. The decay constant 


Gi(q') = i 
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for J - j, and 


amplitude can be written as: 
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(14) 


for J - In Eqs. (13) and (14), Mai. and Fai* are the masses 
and widths of these intermediate N* states, respectively. The 
values used in the present work for Mai* and Tai* are summa¬ 
rized in Table. I. 


TABLE I: Relevant resonant parameters for A* states. The values are 
taken from Particle Data Book [15]. 


A* 

(MeV) 

^ (MeV) 

A(938) 

938 

0 

A(1440) 

1430 

350 

A(1520) 

1515 

115 

A(1535) 

1535 

150 

A(1650) 

1655 

140 

A(1720) 

1720 

250 


On the other hand, we also need to introduce the form fac¬ 
tors for these intermediate off-shell N* {N), which are taken 
as in Refs. [25-28]: 


F(q^) - 


A4 -H (q'^ - M2.)2 ’ 


(15) 


where the cutoff parameter A can be parameterized as 


A - Mat. + /3Aqcd, (16) 

with Aqcd = 220 MeV. The parameter (3 will be determined 
by fitting the experimental data. 

For the background contribution depicted in Fig. 1 (b), we 
construct the amplitude in analogy of Ref. [29]: 

MnoR ^ gNoRv(k2)ey>^u(_kx)-u(_p2)y^y5v{p^)TNoR{s), (17) 

with Tnor{s) — exp(-a( yfs - S/m/)^), where means 

the mass of the final states are summed over. The parameter 
a will be fitted to the experimental measurements, and s is the 
invariant mass square of the system. 

In the phenomenological Fagrangian approaches, the rela¬ 
tive phases between amplitudes from different diagrams are 
not fixed. Generally, we should introduce a relative phase be¬ 
tween different amplitudes as free parameters, and the total 


MNoRe"’”""^ + v{k2)ey^u{ki)-^emllf^ 


-gya+l^ r 


s - + im^Y^ 




i(pN* 


(18) 


where describing the subprocesses i/r(3770) ^ 

are given completely in appendix. 

The differential cross section is given by [30] 
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(Ijlf 2 \Me*e-^ppA^ , , 
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4V(^1-^2)2 


(19) 


and the phase space factor is given by 


^i 03 


1 1 

(2;r)9 8V^ 


|^||p2|^/Q3c/Q2c/nVi 


( 20 ) 


with 2 |A1p averaging over the spins of the initial and 
summing over the polarizations of the final states pp. 

As we can see in the appendix, in the tree-level approxi¬ 
mation, only the products like = gvNN'gnNN' enter in the 
invariant amplitudes. They are determined with the use of MI- 
NUIT, by fitting to the low energy experimental data on mass 
distribution of e^e~ PpA yfs - 3.773 GeV [10]. So 
far we have fifteen unknown parameters: six gisi>, six phase 
angles ^at* and (p^oR, one cutoff f3 in the form factors and two 
parameters gf^oR and a in direct production amplitude Fq. ( 1 7). 
We perform those fifteen-parameter fits to the BFSIII ex¬ 
periment data on the invariant mass distribution at 3.773 GeV 
below 1.8 GeV, and make use of the total cross section infor¬ 
mation in Ref. [10]. Here, we do not consider the invariant 
mass region beyond 1.8 GeV, where contains large contribu¬ 
tion from higher mass N* states and other complicated reso¬ 
nance which decays to pp. In Ref. [9], it was pointed that in 
the case of pp —> AJIip reaction the higher mass N* reso¬ 
nances are needed. Indeed, in the present case, if we go be¬ 
yond 1.8 GeV, we need also the higher mass N* states. On the 
other hand, we did also another calculation including the con¬ 
tributions of higher spin nuclear excited states, A(1675)5/2“ 
and A( 1680)5/2^. It is find that their contributions are quite 
small and the fitted parameters for the other nuclear resonance 
are little changed. Thus, we will not include the contributions 
of this two states in this work. 

We get a minimal = 1.03 with the fitted cut-off pa¬ 

rameter= 6.2 + 3.5. The parameters appearing in direct am¬ 
plitude Fq. (17) are^AfoS = 0.45+0.02, cpNoR - 4.84+0.20 Rad 
and a - 0.84 + 0.02. The other fitted parameters are compiled 
in Table II. The fitted results are shown in Fig. 2 compared 
with the experimental data taken from Ref. [10], where the 
green dashed line stands for the background contribution, the 
orange doted line stands for the nucleon-pole contribution, the 
red line is the full result, and other lines show the contributions 
from different N* resonances. Notice that we have converted 
the experimental event to physical differential cross section 
using the experimental value (Ttotai = 7.71 pb at 3.773 GeV 
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[10]. Our results can describe the two clear peaks around 1.5 
GeV and 1.7 GeV, thanks to the contributions from A^(1520), 
A^(1535) and A^(1650) resonances. The contribution from the 
nucleon pole is small, while the background contribution is 
quite large. 

In Fig. 2, it is interesting to see large interfering effects be¬ 
tween different contributions. At the low Mp„ region around 
1.1 - 1.3 GeV, large cancelation between the nucleon pole and 
the background leads to quiet suppressed spectrum, and the 
bump structure from the nucleon pole just disappears. From 
the two-peak region around 1.4- 1.8 GeV, we can directly 
see that, the background contribution plus N* contribution 
(means without interfering contribution) is not able to reach 
the data peak, it indicates a large enhancement between the 
background contribution and N* contribution thanks to the in¬ 
terfering effect. 


TABLE II: The fitted parameters in the process e*e —> ppif’, 
where = g^( 3770 )wivg;rw. For nucleon, gjv is defined as gjy = 
8ili(iV0)NN8nNN 


N" gA,. (xl0“^) ^A^'Crad) 
A(938) 8.00 ± 0.46 
A(1440) 1.92 ±0.98 6.09 ± 0.38 
A(1520) 0.28 ±0.24 3.74 ±1.07 
A(1535) 1.74 ± 1.34 2.99 ± 0.67 
^■(1650) 1.99 ±0.18 2.17 ±0.19 
^■(1720) 1.14 ±0.63 6.02 ±0.71 



FIG. 2: (color online). The fitted mass spectrum of the process 
—> ppn° at ^fs = ?>.11?! GeV comparing to the experiment data. 
The experiment data are taken from Ref. [10]. The green dashed line 
stands for the background contribution, the orange doted line stands 
for the nucleon-pole contribution, the red line is the full result, and 
other lines show the contributions from different N* resonances. No¬ 
tice that the experimental event is converted to physical differential 
cross section using the experimental total cross section at yj = 3.773 
GeV [10]. 
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FIG. 3; (color online). The typical Feynman diagrams for the process 
pp TT^tf/OnO). 


III. THE y(3770) PRODUCTION IN THE PROCESS 

pp i/f{3no)7T° 


A charmonium plus a light meson n can produced by the 
low energy pp annihilation process. The tree level diagrams 
for the pp tA(3770)7r° reaction are depicted in Fig. 3. It 
is worth to mention that the effect of the N* resonances in 
the cross channel of Fig. 3 has been studied firstly in the 
pp —> 71^ JIlf/ reaction [9]. It was found that the contributions 
from the N* resonances in the pp ^ reaction are im¬ 

portant. In the present work, we extend the model of Ref. [9] 
to the process of the higher charmonium states [i/r(3770) and 
1/^(3686)] production. ' 

The differential cross section of the pp —> reac- 


* We mention that the Regge exchange may be important, unfortunately, the 
information of the Regge propagators are scarce and we hope we can in¬ 
clude the Regge contribution in the future. 


tion at center of mass (c.m.) frame can be expressed as [15] 


do'pp^TP^OnO) ^ 1 1 ^ 3 ™! y jyp 

dcosO 2>2ns \p™\ ^ 


( 21 ) 


where 0 denotes the angle of the outgoing relative to beam 
direction in the c.m. frame, p™ and p™ are the three- 
momentum of the proton and ^/t(3770) in c.m. frame, re¬ 
spectively, while the total invariant scattering amplitude Ai 
is given in appendix using cross symmetry. 

With the parameters determined from the process of 
e^e~ ppTr“, we calculate the total and differential cross 
sections of pp ^ 7r“i/f(3770) reaction. In Fig. 4, we show our 
results for the total cross section of the pp —> 7r°i/r(3770) re¬ 
action as a function of the invariant mass (Ecm) of pp system. 
At Fcm = 5.26 GeV, the total cross section is 0.056 nb, and it 
is under the upper limit of the value obtained in Ref. [10]. 

From Fig. 4, we see that the nucleon pole gives largest con¬ 
tribution, and becomes dominant in the region Ecm > 5.0 
GeV. This is because in the reaction of pp —» 1/^(3770)77*^, 
the four momentum square, of nucleon or other nucleon 
resonance is smaller than 0, and the propagator will in- 
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FIG. 4: (color online). Total cross section of the pp —> n°if/(3n0) 
reaction. The black line is total result, and other lines show the con¬ 
tributions from different N" resonances. 


crease the contribution of nucleon because of its small mass. 
Besides, it is found that the contributions from N* state with 
different quantum numbers have quite different behavior. The 
contributions from A^(1535) and A^(1650) with j de¬ 
crease at Ecm around 4.2 GeV, while the others increase all 
the time. Overall, the total cross section become quiet flat 
while Ecm >4.1 GeV. 

In addition, we also calculate the angular distribution of the 
pp —» n^if/{311Q) reaction at item = 4.0, 4.25, 4.5, 4.75 and 
5.0 GeV. The numerical results are shown in Fig. 5. We can 
see that there emerges an obvious peak at the backward an¬ 
gles (around cos0 ~ -0.8) at Fcm ^ 4.25 GeV produced by 
the contributions of nucleon results in the i/-channel, while 
the larger results at the forward angles is due to the f-channel 
nucleon resonances contributions. 



FIG. 5: (color online). Angular distributions of the pp —> 
reaction with full contribution. 

In Fig. 6, we show the numerical results of the angular dis¬ 
tributions by only considering the contribution from the nu¬ 


cleon pole. We can see that the angular distributions are sym¬ 
metry between the backward and forward angles. Comparing 
Fig. 5 with Fig. 6, we see that, there is a big difference between 
the full contribution and the only nucleon contribution. Our 
model predictions may be tested by the future experiments. 
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FIG. 6: (color online). Angular distributions of the pp —> 
reaction considered only the contribution from the nucleon pole. 


Note that the exchanged nuclear resonances in Fig. 3 are 
far off mass shell, and the form factors for exchanged nuclear 
resonances here should be different with those that have been 
used for the i/'(3770) ^ ppj^^ reaction. We know that 

the form factors can be directly related to the hadron struc¬ 
ture. However, the question of hadron structure is still very 
open, we have to adjust the form factor to fit the experimental 
data, and the hadronic form factors are commonly used phe¬ 
nomenologically [25-28]. The effects of these form factors 
could substantially change the predicted cross sections. Be¬ 
cause of the lack of the available experimental measurements, 
we can not determine the form factors without ambiguities. 
In the present work, we take the same form factors for both 
pp —» t/'(3770)7r*' reaction and e^e~ tl/(371Q) ppjfi re¬ 
action. 


IV. THE IMPLICATION FOR 1^(3686) ^ ppn° AND 
pp —> i/'(3686);r® 

For the process i/'(3686) —» ppn^, we first determine 
the coupling constant g^( 36 & 6 )NN, i-S-, by using the La- 
grangian in Eq. (4), g^(p,6&6)NN can be fitted through the pro¬ 
cess i/r(3686) ^ pp. With the experimental value [15] 
B(i/r(3686) —> pp)- 2.8 x 10“"^, gfpesewA' is determined to 
be 

gip06S6)NN = 9.4 X 10 (22) 

which is consistent with that given in Ref. [4]. 

In Ref. [14], BESIII released the pjr invariant mass spec¬ 
trum of the process ^(^(3686) —> pp^r^ and decay width 
F((/r(3686) ^ pptt") = (1.65 + 0.03 + 0.15) x 10^^ Simi¬ 
lar to the case of i/f(3770), we fit five coupling constants g/^-. 
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FIG. 7: (color online). The fitted pn invariant mass spectrum of the 
process of i/f(3686) —> ppn°- The dashed green curve stands for the 
contribution of the nucleon pole, the solid red line stands for the full 
contributions, and other lines show the contributions from different 
N* resonances. The experiment data are taken from Ref. [14]. 


five phase angles and a cut off parameter /5 to the experimen¬ 
tal data. The fitted results are shown in Fig. 7. Here, one gets 
X^/d.o.f = 2.90 and/? = 3.28 + 2.23, while the fitted coupling 
constants and phase angles are listed in Table III. 

TABLE III: Fitted coupling constants gjv* and phase angles in 
the process i/r(3686) pp7T°, where = g^(ims)NN'gnNN’ 


A* gv.(xl0-^') 

(rad) 

A(1440) 5.10 ±0.86 

3.40 

+ 

0.22 

A(1520) 2.27 ±0.39 

4.96 

+ 

1.10 

A(1535) 0.51 ±0.36 

0.75 

+ 

0.64 

A(1650) 0.76 ±0.19 

5.35 

+ 

0.92 

A(1720) 0.98 ±0.42 

1.77 

+ 

0.99 


In Fig. 7, the dashed curve stands for the contribution of 
the nucleon pole, the solid line stands for the full contribu¬ 
tions, and other lines show the contributions from different 
N* resonances. We see that we can describe the experimental 
data fairly well. Furthermore, we find that the peak between 
1.6 GeV and 1.7 GeV mainly comes from the contribution of 
Ai(1650). 

There also exist quiet obvious interfering effects between 
different A* contributions in Fig. 7. Close to Mp„ = 1.6 GeV, 
comparing the V(1440) contribution to the total contribution, 
one can see the V(1440) contribution is ’’digged out” a valley 
by other N* contributions. In the region of >1.7 GeV, 
the total contribution is smaller than the N{ 1440) contribution, 
i.e., the total contribution is suppressed by interfering terms. 
So, from Fig. 2 and Fig. 7, one can see how important the 
interference effect is. We will not be able to get a good fit 


without interfering terms and arbitrary phase angles. 

Additionally, we also calculated the branch fractions of 
1/^(3686) ^ iN*p + c.c.) from individual intermediate 

N* (or p) state, with the fitted coupling constants listed in the 
Table III. Our results are shown in Table IV. The errors of our 
theoretical results are obtained from the errors of those fitted 
coupling constants of gi^-. We also notice that in Ref. [14] 
BESIII also extracted the corresponding branching fractions 
without considering the interference of different intermediate 
N* (or p) states, which is different from the treatment in the 
present work. Thus, in Table IV we further compare our result 
with the experimental results [14], we see that our results are 
in agreement within errors with that given in Ref. [14]. 

TABLE IV: The calculated branching fractions i^(3686) —> ppn^ if 
considering individual intermediate N" (or N) contribution, and the 
comparison with the experiment values of Ref. [14]. Here, all values 
are in the unit of 10“^. 


Our results The results in Ref. [14] 


N 

7.5 

f\ 49+0.20-F1.78 
^•‘^^-0,20-1.28 

A(1440) 

14 ±4.5 

TJ co+0.25+1.59 
25-0.84 

A(1520) 

2.8 ±0.8 

n ^4+0.054-0.22 
-0.05-0.17 

A(1535) 

2.1 ±3.0 

o '7+0.28+0.99 
^•^'-0.28-0.97 

A(1650) 

4.9 ±2.5 

7 '7^+0.28+1.37 
'^-0.28-1.66 

A(1720) 

1.3 ± 1.0 

1 7Q+0.10+0.24 
^ - '^-0.10-0.71 



FIG. 8: (color online). The cross section of the process pp 
;r'V(3686). The black line is total result, and other lines shows the 
N" contribution. 

With these fitted parameters, we calculate the cross section 
of the process pp —» 7r*^iA(3686) with cross symmetry. The 
results are shown in Fig. 8. One can see that the nucleon 
pole contribution is predominant in the whole energy region, 
while the contributions from other N* states are small. In the 
higher energy region, the nucleon pole contribution is start¬ 
ing to decrease, while the full contribution increases slowly. 





































7 


this behavior resembles the process pp —> 7^1^(3710). Fur¬ 
thermore, it is noticed that, the discrepancy between the total 
result and the nucleon contribution is smaller than the case of 
pp 7rV(3770). 

Finally, we show the angular distributions of the process 
pp —> ;r‘V(3686) in Figs. 9 and 10. Similar to Fig. 5, there 
is a peak in backward angle and a valley close to cos 6 -0. 
Comparing to the angular distribution with the nucleon con¬ 
tribution in Fig. 10, there exits obvious difference, since the 
nucleon contribution only is symmetry while total contribu¬ 
tion is asymmetry. 



FIG. 9: (color online), the angular distribution of the process pp —> 
;r'V(3686). Each line shows the different c.m. energy. 



FIG. 10: (color online), the angular distribution of the process 
pp —> n°tf/{36&6) only considering the nucleon contribution. Each 
line shows the different c.m. energy. 


V. DISCUSSION AND CONCLUSION 

We have studied the PP^r*’ at 3.773 GeV c.m. en¬ 

ergy and pp —> 71^1/^(3770) reaction within an effective la- 


grangian approach. The —> ppn^ process is a good plat¬ 

form to study excited N* nucleon resonances. We consider 
contributions from nucleon pole and five well established N* 
states. First, we perform a;^f^-fit to the experimental data on 
the mass distribution of the —> pp7r“ , from where we 

obtain the couplings of i/^(3770) to these N* states. It is shown 
that we can describe the experimental data quite well. In par¬ 
ticular, the two bumps around 1.5 and 1.7 GeV can be well 
reproduced. We also find that the contribution of the nucleon 
pole is small comparing to the background contribution, and 
there exists large cancellation in low Mp„ region. 

Second, based on our results of the ^ we 

study the pp —> 7r“i/^(3770) reaction with cross symmetry. We 
evaluate the total and differential cross sections of the pp 
7r“i/7(3770) reaction. The nucleon pole gives largest contri¬ 
bution to the pp —> 7r°i/^(3770) reaction close to threshold. 
However, the interference terms between nucleon pole and the 
other nucleon resonance affects significantly and could change 
the angle distributions clearly. Our studies provide valuable 
information to future experimental exploration the i/r(3770)7r“ 
production through the pp interaction. 

Additionally, we also study the tl/(36^6) production through 
the process pp —> ;7r‘'i/^(3686). Similarly to the case of 
—> 1/^(3770) —> pp7r“, we study firstly the decay pro¬ 
cess of 1/7(3686) ^ ppn^ to extract the parameters we needed. 
Then we study the pp 7r*'tf7(3686) reaction. We find that 
the contribution from the nucleon pole is dominant, while the 
angular distributions show a quite discrepancy induced by the 
N* states. 

We hope and expect that future experiments at PANDA will 
provide a test to our model and give more constraints on our 
theoretical study. 
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Appendix: Scattering amplitudes of the subprocess 

1/7(3770) —> pp7T° and the process pp —> Tfiij/OinO). 

The tree level diagrams of the subprocess iJ/{377Q) —> ppTr'^ 
is depicted in part of Pig. 1 (a). According to the feynmann 
diagrams shown in Pig. 1, the scattering amplitudes M.jp with 
a exchanged (including N) are given by: 
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M 


gjiNPi 


-gVNPn^{Pl)u{p2) 


- 2mi\i 

^1- - g7iNSngVNSi,^(Pl)l<(P2) 


-pt + niM’ Plu + 

+ nim)- — —r^Fiu) 


t — m 


757p- 


N’ 

fft + mN> 


V(P3), 


Ml* = ^^^^igVNPu^{P\XP2) 
2 niM 

Ml- = ^^^^^gVNDn^iP\)u{P2) 


75- 


-0, + mn- 


F{t) + 


Plu + »%• 


757pF(u) 


V(P3), 


t — m 


Plu + m«. 

2 G,,{-p,){ip\)F{t) + - —Gy^ipuIrsFiu) 


w 

-[f, + niN' 


t — m 


tiu + niMt 

2 G,y(-p,)j5(i[l4)iipl)F(t) + 75(iff4){ipl)- - —G,,{pu)F{u) 

w- “ - 


(23) 

(24) 

V{P3), (25) 

v{P3) (26) 


with Pi ^(pi- p 2 r and u ^ pi ^ (pi - P 3 ) , and G^y is 

1 


G,y(p) - ( g,y+ ^y.Jy . 


1 , 

i7pPv-7yPp)-^ n—^)- 


(27) 


For in Eq. (18), just drop the polarization vector e^{p\). 

For pp TiP^{'illO) reaction, the scattering amplitudes can be easily obtained just applying the substitution to Eqs. (23-26): 

El ^ -P3, P2 -» -P2, P3 -Pu Pt -» -Pu- (28) 

The amplitude of the process i/'(3686) ^ pp^^ and pp —» 7r**iA( 3686) is exactly the same. _ 
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